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Abstract: Photons couple to Axion-Like Particles (ALPs) or more generally to any pseudo 
Nambu-Goldstone boson in the presence of an external electromagnetic field. Mixing be- o 
tween photons and ALPs in the strong magnetic field of a Gamma-Ray Burst (GRB) jet 
during the prompt emission phase can leave observable imprints on the gamma-ray polar- 
ization and spectrum. Mixing in the intergalactic medium is not expected to modify these 
signatures for ALP mass > 10~ 14 eV and/or for < nG magnetic field. We show that the 
depletion of photons due to conversion to ALPs changes the linear degree of polarization 
from the values predicted by the synchrotron model of gamma ray emission. We also show o 
that when the magnetic field orientation in the propagation region is perpendicular to 
the field orientation in the production region, the observed synchrotron spectrum becomes 
steeper than the theoretical prediction and as detected in a sizable fraction of GRB sample. 
Detection of the correlated polarization and spectral signatures from these steep-spectrum 
GRBs by gamma-ray polarimeters can be a very powerful probe to discover ALPs. Mea- 
surement of gamma-ray polarization from GRBs in general, with high statistics, can also 
be useful to search for ALPs. 
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1. Introduction 

Axions represent the most convincing and elegant solution to the strong CP problem Q. 
They are the pseudo Nambu-Goldstone bosons of a global U(1)pq symmetry j|, [|. Axion- 
Like Particles (ALPs) can be understood as generalizations of the axions, and appear 
generally in theories beyond the Standard Model (SM) of Particle Physics. Axions and 
ALPs couple to photons in the presence of an external electromagnetic field j|] . For axions 
the strength of the coupling g ai is inversely proportional to the energy scale M at which 
the U{\)pq symmetry is spontaneously broken, and is directly related to the particle mass 
m a . For ALPs there is no general relation between the energy scale of the new physics 
beyond the SM and the ALP mass, and therefore in the following we shall consider the 
parameters g ai and m a to be independent of each other. ALPs may be copiously produced 
in the early universe, either thermally Q or non-thermally @, providing a possible (sub) 
dominant (hot) dark matter candidate. 

Mixing of photons and ALPs in an electromagnetic field results in photon- ALP con- 
version and a change in photon polarization states. The former effect has been extensively 
exploited to search for ALPs that are created in the Sun, travel to the Earth as ALPs 
and convert to ~ keV photons in the magnetic field of a laboratory experiment. From 
non-detection of such photons, the CAST experiment has reported a lower bound on the 
ALP energy scale of M > 1.1 ■ 10 10 GeV, which translates to a constraint on the photon- 
ALP coupling of g ai < 8.8 • KT 11 GeV" 1 for ALP masses of m a ~ 0.02 eV §. These 
constraints exclude a region in the g^-ma parameter space. The same conversion mech- 
anism is used by the ADMX experiment to search for ALP dark matter that converts to 
microwave photons [p. 

Observation of supernovae (SNe) la dimming has been suggested as a possible signa- 
ture of photon-to-ALP conversion, thus depleting the photon flux, in the Inter-Galactic 



Magnetic Field (IGMF) g, [l(| 0, ||, |T|, @. This is an alternative to the standard inter- 
pretation by a dark energy fluid that is responsible for recent accelerated expansion of the 
Universe, making the distances of the SNe la larger. Search for circular polarization in the 
Cosmic Microwave Background (CMB) data has been proposed as probe of photon-ALP 
mixing in the IGMF [15]. Mixing in the IGMF has also been considered as a possible 



mechanism to produce ultra-high energy cosmic-ray events, assumed to be photons which 
are not attenuated while in their ALP states and while propagating from distant sources to 



the Earth [16|. A similar mechanism has been proposed to search for photon-ALP conver- 



sion effects in the GeV-TeV 7-ray fluxes from distant active galactic nuclei [[17], |1S|, [L9|, g0|. 
Detection of these fluxes at very high energies may provide hints of photon-ALP mixing, 
which would be absorbed by the Extragalactic Background Light (EBL) otherwise (see e.g. 
Ref. |H). 

Here we study the observational consequences of a photon-ALP coupling on GRB 
photon polarization and fluxes at low energies, in the ~ keV-MeV range, arising from high 
magnetic field in the GRB jet. 1 A search in this energy range has several advantages: 
(i) GRBs, the most powerful explosions in the Universe, release upwards of 10 53 erg of 
isotropic-equivalent 7-ray energy, mostly in the ~ keV-MeV range pq] . Thus we have the 
most powerful photon beam at our disposal to investigate the effect, (ii) Unlike TeV 7 rays, 
MeV photons are not attenuated in the EBL and detection of any effect due to photon- 



ALP coupling does not depend on the EBL models (see e.g. Refs. ]26|, |27j, [28|, [29|, pC|1 ). 
(iii) Since photons are converted to ALPs in the high magnetic field of the GRB jet, 
ALPs may not convert back to photons while propagating in the IGMF and in the galactic 
magnetic fields 2 . Search for photon-ALP mixing in the photon polarization data has been 
suggested for GRBs in the past, both in the strong magnetic field of the GRB |3l| and 
in the IGMF ||2| (see also Ref. |33|]). We have studied effects of photon-ALP mixing on 
GRB 7-ray polarization using a realistic emission model, namely synchrotron radiation by 
relativistic electrons in the strong magnetic field, either advected from the GRB central 



engine |34|, |35| or generated in the shocks |36|, 37, j^, |3f| or both. 

In the synchrotron model, which is also the leading model for observed ~ keV-MeV 
7-ray emission, a population of electrons are assumed to be injected as a power-law above a 
minimum particle Lorentz factor in the magnetized plasma with an optical depth less than 
unity. The peak of the observed energy spectrum (E 2 dN/dE), typically in the ~ 0.1-1 MeV 
range, is identified with the characteristic synchrotron frequency from the electrons with 
the minimum Lorentz factor in the comoving GRB jet frame, boosted by the bulk Lorentz 
factor of the jet. Synchrotron radiation is partially polarized with a linear polarization 
degree of ~ 50% at frequencies much lower than the characteristic frequency, reaching 
~ 70% at the maximum |40[| . We model the initial polarization states of the observed 
photons in the ~ keV-MeV 7 rays as from the electrons with minimum Lorentz factor, 
according to synchrotron radiation theory in the comoving frame. The effect of photon- 



1 GRB pseudo-Goldstone boson emission and its subsequent conversion to electromagnetic energy was 
proposed as a possible mechanism for the observed GRBs p2| , [23] , pi| . We do not address such a possibility 
in the present study. 

2 They can, however, convert back to photons in a suitable laboratory experiment. 



ALP mixing then changes the observed polarization from the expected pattern. 

To date, prompt 7 ray polarization has been measured from only a handful of GRBs, 
most notably an (80 ± 20)% linear polarization from GRB 021206 by RHESSI @. How- 
ever, these measurements are statistically inconclusive and suffer from large systematic 
uncertainties (see e.g. Ref. p2| ). Gamma- Ray Burst Polarimeter (GAP), sensitive in the 
50-300 keV range, aboard the recently launched IKAROS Solar Sail is one of the new gen- 
eration of instruments to measure 7 ray polarization [43]. A number of satellite missions 
such as the Advanced Compton Telescope (ACT) [44], Gamma-ray Burst Investigation via 
Polarimetry and Spectroscopy (GRIPS) |f45| , and Polarimeters for Energetic Transients 
(POET) E(| are also being planned to measure 7 ray polarization in the keV-MeV range. 
These experiments are expected to measure GRB polarization with a high statistical signif- 
icance and have been shown to be excellent tools to test the synchrotron emission models 
(see e.g. Ref. p?J). Eventually, these broadband polarimeters will be able to detect de- 
viations from the standard synchrotron polarization pattern. Such frequency-dependent 
deviations in the polarization pattern could be explained in terms of photon- ALP mixing. 

The polarization pattern induced by photon- ALP mixing can be accompanied with a 
detectable change in the 7-ray spectral slope, due to a depletion of preferentially low energy 
photons that convert to ALPs in the GRB jet. Indeed a specific prediction of the GRB 
synchrotron model is that, below the peak energy the spectrum can not be harder than 
the photon index 07 = —2/3, where dN/alE oc E a ~> , a limit that arises from synchrotron 
theory of radiation from a single particle ]^0|, ^] . Observed variation of the GRB low-energy 
spectra softer than this limit may be explained as cooling effect on the electron spectrum, 
producing a 7-ray spectrum as soft as dN/dE oc E~ 3 / 2 (see, e.g., Ref. [§9|). Majority of 
bright GRBs, detected by the Burst And Transient Source Experiment (BATSE) aboard 
the Compton Gamma Ray Observatory, for which good spectral data are available p0[ 



falls within the synchrotron limit of —2/3 > 07 > —3/2. However a significant (~ 20%j 



fraction violates the "synchrotron death line" of 07 = —2/3 [51], and a "harder when 
brighter" tendency is present in the data. The same effect has been detected in time- 
integrated and time-resolved spectra from joint observations by the Burst Alert Telescope 
(BAT) aboard Swift and by the Wide band All-sky Monitor (WAM) aboard Suzaku J52|, 
and most recently by the Gamma-ray Burst Monitor (GBM) aboard the Fermi Gamma-ray 
Space Telescope |33|. We predict that polarization measurements of these steep-spectrum 
GRBs can shed light, or even lead to discovery of ALPs. 

The structure of the paper is as follows. We review the photon- ALP mixing phenomena 
in Sec. 2 and apply this formalism to the GRB jet and synchrotron emission model in Sec. 
3. We discuss our results in Sec. 4 and conclude our study in Sec. 5. 



2. Photon-ALP mixing and conversion probabilities 



We follow here the photon-axion/ALP interaction formalism from Ref. || (see also Ref. [p 7 ^]). 



The lagrangian for the photon- ALP system is given by 3 
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where F^ v is the electromagnetic field tensor, F^ u = \e^ upa F pu is its dual, a is the fine- 
structure constant and m e is the electron mass. The second term in Eq. ([O]) is the 
Euler-Heisenberg effective Lagrangian, which accounts for one-loop corrections to the clas- 
sical electrodynamics. The third and fourth terms in Eq. ( |2.1| ) are the Lagrangian terms 
describing the ALP field a with a mass m a . The last term is the photon- ALP interaction 
lagrangian, which, in terms of the external electromagnetic field, reads 
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(2.2) 



Here <? a7 is the photon- ALP coupling constant, E and B are the electric and magnetic 
fields respectively. 

The evolution equations for a mono-energetic photon/ALP beam with energy oo prop- 
agating along the z direction in an external and homogeneous magnetic field transverse 
(B7 1 ) to the beam direction (i.e. in the x-y plane) are given by: 
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Here A± and An are the two photon polarization components (both in the x-y plane) 
perpendicular and parallel to the external magnetic field B^, respectively. The plasma 
term in the equations of motion arises due to the presence of electrons in the media, giving 
an effective mass to the photons, and is proportional to the electron number density n e . 
The critical magnetic field is defined as -B cr i t = m 2 /e = 4.414 • 10 13 G, where e is the 
electron charge. In the limit where ui ^> m a , the evolution of the system can be linearized 
in the form of a first order differential equation 4 

Here M. is a mixing matrix of the axion field with the photon polarization components, 
and is given by 
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We adopt the natural unit convention h = c = 1. 
4 We follow the notation adopted in Ref. ]32]]. 
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The elements of M can be expressed as A_j_ = 2Aqed + A p i, An = (7/2)Aqed + A 



pi 



following Ref. [32], and we provide their reference values relevant in our case below 



A = \ g Bt ~ 1.32 . 10- 11 ( T ) ( J%-) cm- 1 , 

7 2 ya < \ 8.8 ■ lO-HGeV -1 / V 106(3 / 

A ^_^_ 2 . 53 . 1(r i3(_) (_^L_) ( 2 . 6) 



The plasma frequency is defined as u p \ = \J ' Airan e /m e = 3.71 • 10 14 y / n e /cm 3 keV. 
Notice from Eqs. (2.4) and (|2.5|) that the component of the photon beam polarization 



perpendicular to the By field, A±, will decouple from the evolution of the photon- ALP 
system. In other words the ALP couples only to the An polarization component. 

A generalization of the scenario discussed so far is when Bt makes an angle £, < 
£ < 2tt, with the y axis in a fixed coordinate system. A rotation of the mixing matrix 
[Eq. Q2.5D 1 in the x-y plane then leads to a new form and the evolution equation of the 
photon- ALP system reads 




( A^cos 2 £ + Ay sin 2 £ cos£sin£(A|| - A_l) A a7 sin£ 
cos£sin£(A|| — Aj_) Aj_ sin 2 £ + An cos 2 £ A a7 cos£ 
^ A a7 sin£ A a7 cos£ A a 



(2.7) 

If there are more than one magnetic field domain present in the problem, then Eq. ( [2.7] ) 
needs to be solved for each domain with appropriate initial conditions. Under the assump- 
tions that all the domains in a particular environment (constant n e and the same initial 
conditions for the fields at z = 0) have identical coherence lengths and magnetic field 
strengths, and only the orientation of the magnetic field Bt in each domain is random, 
then the average effect can be calculated by randomly varying £. We mainly consider the 
scenario where photons are created at z = 0, at source, and cross Bt field domains where 
(i) £ = or 7r/2 in each domain, and (ii) £ is random. In both cases each photon crosses 
only one coherence length width in the z direction. The final beam consists of contribu- 
tions from all domains. This is different from propagation of the beam in the intergalactic 
medium where each photon/ ALP crosses many IGMF domains and the initial conditions 
change each time the beam enters a domain (see, e.g. Ref. [p2[ ). 

In analogy with two-family neutrino mixing, the conversion probability of An into 
ALPs after traveling a coherence length L and for £ = reads 

P a7 =sin 2 20sin 2 f^^y (2.8) 



where the oscillation wave number is A osc = w (A a — Ay) 2 + 4A 2 7 and the mixing angle is 
9 = (1/2) arctan[2A a7 /(Aii — A a )]. From Eq. fl2.8| ) it is possible to infer the energy range 



in which the conversion probabilities are approximately energy independent and mixing 
effects will be maximal {0 ~ vr/4) for ujl < uj < cjff, where the low and high critical 
energies, respectively, are given by [19, 32] 
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We focus on the specific problem of photon- ALP mixing in the GRB jet and the impact of 
the photon- ALP conversions in the observed photon spectrum in the next section. 

3. Gamma-ray emission and conversion to ALPs in the GRB jet 

Synchrotron radiation from relativistic electrons that are accelerated in the GRB jet, either 
due to dissipation of the jet kinetic energy (internal shocks of plasma shells) p6[ or mag- 



netic flux [34] from a central engine, is believed to be the dominant mechanism to produce 
observed 7 rays in the keV-MeV range. In the internal shocks model the conversion of jet ki- 
netic energy to 7 rays takes place at a radius R 2T 2 ct v ~ 2.7-10 13 (T/300) 2 (i„/10 -2 s) cm, 
which can vary widely depending on the jet bulk Lorentz factor V and the 7 ray flux vari- 
ability time scale t v . The jet kinetic energy is typically estimated from the observed 
isotropic-equivalent 7-ray luminosity L 7 and assuming that a fraction e e of the kinetic en- 
ergy is converted to relativistic electrons which promptly radiate most of their energy to 7 
rays. Random magnetic field in the GRB jet is believed to arise when a fraction e# of the jet 
kinetic energy is converted to the magnetic field energy in the shocks (see e.g. Rcf. p9| , Q). 
An average value of the magnetic field and electron density can be estimated, in the jet 
comoving frame, as B ~ 5 • 10 4 (e B /e e ) 1 / 2 (L 7 /10 52 erg s- 1 ) 1 / 2 (r/300)~ 3 (t v /10- 2 s)" 1 G 
and n e ~ 2 • 10 8 e^^/lO 52 erg s' 1 )^ /300)- 6 (t v /10~ 2 s)" 2 cm" 3 (see e.g. Ref. §5|). 

Strong magnetic field from the central engine can also be present in the GRB jet. The 
toroidal component of the magnetic field of a magnetar with surface magnetic field Bo at 
#0 ~ 10 6 cm drops to a value B = B {R /R) ps 10 8 {B /W 15 G)( j R/10 13 cm)" 1 G at a 
dissipation radius R. The magnetic field from the central engine is globally ordered in 
the emission region. The coherence length scale of the random magnetic field can be as 



small as the plasma skin depth [39, 54], however efficient conversion of the shock energy 
to 7 rays requires a length scale of the order of the comoving width of the plasma shell 
(Ai?) ps Tct v ~ 9 • 10 10 (r/300)(^/l(r 2 s) cm. Because of relativistic beaming, only an 
angular size scale 1/r of the jet surface is viewable. Note that, this also corresponds to 
a maximum length scale (Ai?) over which the random magnetic field can be fully ordered 



due to causality [54]. The jet half-opening angle 0j et is much larger than 1/r during the 
prompt 7-ray emitting phase. Both the ordered and random magnetic fields are mostly 
perpendicular to the jet axis, which is assumed along the z direction. 

Synchrotron radiation from the visible patch of the jet surface can reach the maximum 
polarization degree, w 50%-70%, if the magnetic field is fully ordered in the patch and 
r#jet 3> 1- Intrinsic curvature of the field, for example in case of toroidal field configuration, 



in a large visible patch can reduce the maximum polarization degree to ~ 40% |_, 47]. 



Smaller scale random magnetic field, if dominant, can also reduce the net polarization 



degree [54, 47]. We explore both the ordered and random field scenarios to calculate 
photon- ALP mixing in the GRB jet. Moreover, the emission region and propagation region 
of the photons can be separated with different magnetic field strengths and orientations 
(i.e. £ 0). Faraday rotation of the polarization plane can be important for synchrotron 
radiation _|] only for a substantial magnetic field component parallel to the beam direction 
(along the z axis) and below the optical frequencies, both situations are outside the scope 
of this paper. Mixing of the An and A± components in our scenario takes place through the 



off-diagonal terms in the mixing matrix [Eq. ( |2.7| )1, due to £. Additional ordered magnetic 



field (e.g. in the wind of the progenitor star) surrounding the GRB jet [56], if present 
and is sufficiently strong, can modify some of the polarization effect that we explore here. 
However we ignore that for simplicity. 

The two photon polarization components in synchrotron radiation can be written in 
terms of the Bessel functions (see e.g. Ref. |Q) as 5 



AM 



from a single electron with Lorentz factor 7 e gyrating in the B field. Here 6 e is the angle 
between the line of sight and the plane containing the electron trajectory. The characteristic 
synchrotron frequency, in case 9 e — > 0, is given by 

3 B sin r] n 

Wc= o^ le m e, (3.2) 

^ -Dcrit 

where rj is the pitch angle between the electron's velocity and B. The intensity of syn- 
chrotron radiation is given by 

^L = ^ { \ M ^ + M ), (3.3) 

and the emitted radiation is concentrated in a solid angle dfl = 2tt sinr/ d9 e . The power 
emitted per unit frequency is calculated by dividing the intensity with the orbital period 
of the charge, T = 2n~i e m e /eB, after integrating over the solid angle as 



The degree of linear polarization for a mono-energetic electron is given by [4C] 



_ Pj.M-.PliM 
Ul = PxM + AM ' (3 ' 5) 



5 Note that Ref. Km uses exactly the opposite convention for the polarization components. 



where P±(uj) and P\\{w) are the powers emitted per unit frequency in directions parallel 
and perpendicular to the magnetic field, and can be calculated from Eq. (j3~|). 

The total synchrotron power from a distribution of electrons 6 can be calculated by 
performing the convolution of the power from each electron and by integrating over 7 e . In 
the keV-MeV range of our interest, however, 7 rays from GRBs are modeled as synchrotron 
radiation from the shock-accelerated electrons of a minimum Lorentz factor 7e im . The 
observed peak photon energy in the E 2 {dN/dE) energy spectrum (often denoted as EF{E) 
or vF u ) corresponds to the characteristic photon energy in Eq. ( |3.2|) , after multiplying by 
aT/(l + z) factor, as E pk ~ 3.5 (1 + z)~ x (B sin^/lO 6 G)(7 eim /10 3 ) 2 r 30 o MeV. The typical 
GRB redshift is z ~ 1-2. Higher energy photons, but not too far above -Epk, can be 
modeled as synchrotron radiation from a power-law distribution of electrons above 7 ejm 
and do not couple to ALPs in our present study. 

To explore photon- ALP mixing in the GRB jet environment, we solve the field evolution 
equation [Eq. (|2.7D 1 with mixing matrix elements [Eq. ( |2.6[ )1 derived from GRB environment 
parameters, and with initial electromagnetic field input from Eq. ( |3.1[ ). Note that the 
comoving frame values for the GRB parameters are used to evaluate photon- ALP mixing, 
and the resulting effect show up in the comoving frame frequency to. The observed photon 
energy is E = uT/(l + z). We calculate the effect of photon- ALP mixing on the polarization 
pattern by using A± and An from solutions of the evolution equation [Eq. (2.7)] to find the 
linear degree of polarization as 

„ _ P±,alp(w) - P||,alpH ,„ r , 

^J_,ALpM +P||,ALPM 



and compare with Eq. ( |3.5|) , without photon- ALP mixing. We also define a flux modifica- 
tion factor, from Eq. ( |3.4[) , as 

p = P(co) ALP /P(u) , (3.7) 

which shows any deviation from the synchrotron spectra due to photon- ALP mixing in the 
GRB jet. We discuss results from our investigation next. 



4. Results and Discussion 

For the nominal values of the GRB parameters B>t = 10 6 G, n e = 10 8 cm -3 , L = 10 11 cm, 
and for the photon-ALP coupling constant g ai = 8.8 • 1CP 11 GeV" 1 which is very close 
to the current CAST limit 0; strong mixing of photons and ALPs takes place in the 
GRB jet when m a < ^2g ai ujB T < lCT^uVkeV eV, from the condition A 2 < 4A^ 7 . 
Indeed the photon-ALP mixing term A a7 dominates other terms [Eq. Q2.6|) ] in the mixing 
matrix for the nominal GRB parameters, and A osc ~ 2A a7 ~ L~ l [Eq. ( |2.8[ )]. The mixing 
angle 9 is also maximized in this (Ay — A ) < A a7 . The off-diagonal rotation term 

oc (Ay- Aj_) = (3/2)A QED ~ 2- 10 _12 (w/keV) cm" 1 is small at low oj for the nominal GRB 
parameters, but can become significant at high to. Thus it is important to keep all terms in 



See e.g. Ref. Kft for power-law distribution of electron Lorentz factor. 



the mixing matrix and solve the evolution equation [Eq. (|2.3| )1 numerically with frequency- 
dependent initial conditions from Eq. ( |3,1[ ). Photon- ALP conversion mostly takes place in 
a broad observed energy range of E ~ (12-560)(r/100)(l + z) -1 keV [Eq. ( |2.9| )1 for our 
reference parameters. 




Figure 1: Linear photon polarization with and without ALP mixing in the GRB jet for the nominal 
GRB parameters Bt = 10 6 G, n e — 10 8 cm -3 and L — 10 11 cm. We used a photon-ALP coupling 
parameter value g ai = 8.8 • 10~ n GeV -1 along with ALP mass m a — 10" 7 eV. The GRB is assumed 
to be at redshift z = 2 with a jet bulk Lorentz factor L = 100. Also the synchrotron emission from 
the GRB is assumed to peak at « 660 keV in the observer's frame. The polarization degree without 
photon-ALP mixing is shown as the black dashed line obtained by solving the evolution equation 
[Eq. (2.7)]. The solid gray line is the expected polarization from synchrotron theory. The role of 
the final A\\ and A± are interchanged from the initial configuration while £ changes from (blue 
dot-dashed line) to tt/2 (red dotted line). Total polarization from many identical domains but with 
random £ is also shown (purple dashed line). 



Figure |l] shows the effects of photon-ALP mixing in the GRB jet with nominal param- 
eters as mentioned above with T = 100 and z = 2. The peak of the synchrotron radiation 
is assumed at uj c = 20 keV in the comoving GRB jet frame (2 MeV in the rest frame of the 
source or ~ 660 keV in the observer's frame). The initial polarization obtained by numer- 
ically solving the evolution equation [Eq. fl2.7|) l, without photon-ALP mixing, is plotted 
with the black dashed line, which agrees with theoretical expectation (solid gray line). The 
results for photon-ALP mixing are plotted for two cases, £ = (blue dot-dashed line) and 
£ = tt/2 (red dotted line). The change in polarization from the £ = case to the £ = ir/2 
case can be understood as the magnetic field orientation in the initial production region 
and propagation region being aligned parallel with each other in the former case and being 
aligned perpendicular to each other in the latter case. In other words, as an inspection of 
the mixing matrix in Eq. (|2.7| ) reveals, the Ah and A± in the final states are interchanged 
from the initial configuration for £ = tt/2. Observations in limited energy bands, however, 
can not distinguish between the two extreme cases and is expected to be intermediate, since 



polarimeters measure the absolute degree of polarization, highly correlated with £ within 
a unique energy band. On the other hand, a change in polarization degree in different 
energy bands, different from the synchrotron radiation pattern, can be used to search for 
photon- ALP mixing signature. 

Time-resolved measurements over small intervals and around the pulses in the GRB 
light curves are important to ensure that emission from only a small bright spot, in which 
we assume the magnetic field to be fully coherent, of the jet surface contributes in each 
case. For longer exposure, contributions from many domains (assumed identical) on the 
jet surface can contribute. This case, where each domain is assumed to have completely 
ordered field within and only the orientation of the magnetic field direction £ is assumed 
random, is also shown in Fig. [l] with the purple dashed line. 
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Figure 2: Phenomenological GRB energy spectrum with and without photon- ALP mixing with the 
same parameters used in Fig. |l| The spectrum for no photon- ALP mixing (solid thick black line) 
is plotted using the Band spectrum with peak photon energy -Epk.Band = 500 keV, and low (high) 
energy power-law slope aeand ~ —0.6 (/?Band = —5/2). Asymptotically the low energy power-law 
slope coincides with the expected spectrum with a 7 = —2/3 from synchrotron theory (thin solid 
gray line). The effects of photon- ALP mixing arc plotted by multiplying the Band spectrum with 
the suppression factor in Eq. (3/7) for the £ = (blue dot-dashed line), £ = 7r/2 (red dotted line) 
and random £ (purple dashed line) cases. As can be seen, the observed spectra for the £ = 7r/2 case 
can be steeper than the synchrotron spectrum in a limited energy range (thin dashed gray line). 
Also plotted are the energy bands in which Swift BAT and Fermi GBM instruments are sensitive. 
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Figure shows the effects of photon- ALP mixing on the GRB spectrum. A phenomeno- 



logical photon spectrum (dN/dE), called the Band spectrum 57], is plotted (thick black 
curve) with the low-energy index aeand ~ —0.6, high energy power-law index /Seand = —5/2 
and a peak energy -E p k,Band = 500 keV. The GRB is assumed to be at z = 2 as in Fig. || 
with all other parameters for photon-ALP mixing the same as those used for Fig. |T]. We 
assume that synchrotron radiation from minimum energy electrons dominates below the 
start of the high energy power-law part of the spectrum at (2 + aBand)-^pk,Band ~ 700 keV, 
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similar to oj c in the jet comoving frame. The effective low energy power-law index |pl| , 
corresponding to the synchrotron theory, is then a 7 ~ —2/3 (i? 4//3 in the E 2 dN/dE or vF v 
spectrum) as plotted with a thin solid gray line. Approximately 20% of the GRB spectra 
are steeper than this "synchrotron death line" , and in ~ 5% of the time-resolved spectra 



of bright GRBs the spectral deviation is statistically significant [ 50 1 . As shown in Fig. 



photon- ALP mixing for £ = ir/2 case (red dotted line) can change the low-energy spectrum 
to as steep as dN/dE oc E° (thin gray dashed line) from the synchrotron model, depending 
on the parameters we used. The change in the spectrum is not as significant, however, for 
the £ = and random £ cases. 

The £ = 7r/2 case should be less frequent in nature as evidenced by the fraction of GRB 
spectra that violates the "synchrotron death line". High polarization degree, up to 100%, 
is expected in these cases (Fig. ||). Indeed the peak-resolved spectra of GRB 021206 with 
(80 ± 20)% polarization Q show low-energy index as hard as «Band = —0.42 ± 0.05 fl58fl . 
A larger sample of GRBs with correlated high polarization and steep low-energy spectrum 
detected with future polarimeters will be instrumental to probe the photon-ALP mixing 
in GRB jets. Other explanation of steep spectrum by black body, jitter radiation, inverse 
Compton scattering etc. (see e.g. Ref. [j59|) do not generally change the polarization pattern 
the way photon-ALP mixing does and as we discussed here. 
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Figure 3: Contour plots of the flux suppression factor in the ALP parameter plane with fixed GRB 
parameters (left panel), and in the mixed ALP-GRB parameter plane (right panel). The outer, 



middle and inner contours depict the regions in which the flux suppression factor, see Eq. (3.7), 
reaches 70%, 50% and 40%, respectively. For m a < 10 -6 eV, the effect becomes independent of the 
ALP mass and is restricted to a rather narrow range of g a ^B values. 



Figure y, left panel, depicts the flux suppression factor p, see Eq. (|3.7|) , in the m a -g ai 



plane with fixed GRB parameters (Bt = 10 6 G, 



n, 



10 8 



cm 



and L = 10 cm). The 



outer, middle and inner contours depict the regions in which the flux suppression factor 
reaches 70%, 50% and 40%, respectively. 7 The suppression effect has been averaged in the 
7 A smaller suppression factor corresponds to a larger conversion probability of photons to ALPs. 



(10-500) keV energy window. A flux suppression factor < 40%, requires a photon ALP 
coupling parameter close to the current CAST limit. As previously discussed, the mixing 
angle is maximized when the photon- ALP mixing term dominates the system evolution, 
and, for sufficiently small ALP masses (m a < 10~ 6 eV), the effect becomes independent of 
the ALP mass. In the pure axion cold dark matter (CDM) scenario, if the PQ symmetry is 
restored after inflation, a lower mass bound m a > 10 -5 eV applies in order to not overdose 
the universe (see Ref. [6C] and references therein). However, if inflation takes place after 
the PQ transition, much smaller values for the CDM axion mass are still allowed, see 



Refs. j6l], |62|, |63|, |64|. The ALP case is more complicated, since these particles constitute 



a dark matter candidate only under certain conditions. For instance, if ALPs couple 
exclusively to photons, they are excluded as CDM candidates |||. Consequently, no lower 
ALP mass bound is shown in Fig. ||| (left panel), since the role of ALPs as CDM particles 
depends highly on the underlying theoretical model. 

Figure ||, right panel, depicts the flux suppression factor p same as in the left panel. 
The contours are plotted in the g a ^B and \m? a — oj^] plane, where both the quantities are 
closely related to the wave numbers [Eq. ( |2.6| )]. Comparing the ranges of g ai and g a ^B 
values from the plots, significant (< 70%) flux suppression takes place for B ~ 4 • 10 5 - 
3 • 10 6 G in the observed keV-MeV range. Thus detection of photon-ALP mixing effect 
in GRB data can, in principle, be used to probe the magnetic field value in the GRB jet, 
which is somewhat uncertain. Note that a much higher, ~ 10 9 G, field with a ~ 10 6 cm 



coherence length, corresponding to the neutron star radius used in Refs. [31, ^] gives no 



photon-ALP mixing effect in the GRB jet. However, the radius of 7-ray emission region 
is likely to be large to avoid e + e~ pair creations by the photons and thermalization. The 
magnetic field in the jet is thus likely to be small, typical to the values that we used, at 
this large radius. 

Mixing of photons with ALPs, for propagation in the IGMF with generally assumed 
magnetic field -Bigmf = 1 n G and particle density n e = 10~ 7 cm -3 , takes place in 
the frequency range [see Eq. fl2.9| )] ujl ~ 10 7 (m a /10~ 7 eV) 2 (i?iGMF/nG) _1 GeV and 
ujh ~ 6 • 10 9 (i?iGMF/nG)~ 1 GeV for the same g ai parameter from the CAST limit. 
The contribution of plasma frequency to Wi becomes dominant for ALP mass m a <C 
10~ 14 y^e/lO -7 cm -3 eV from the condition uj p \ 3> m a in Eq. (2J3). The correspond- 
ing ujl ~ 0.1(n e /10~ 7 cm~ 3 )(i?iGMF/nG) -1 keV becomes constant. The oscillation wave 
number is A osc « A a7 in this asymptotic range, and the oscillation probability [Eq. (g~J)] 
is P ai « (A Q7 L) 2 « 2 • 10" 3 (#iGMF/nG) 2 (L/Mpc) 2 for Mpc scale coherence length. Thus 
photon-ALP mixing in the IGMF can be important over Gpc scale source distance and 
wash-out the source signature only if the IGMF is of the order of nG and the ALP mass 
is smaller than 10 -14 eV. This result is compatible with mixing effect in the IGMF for 
ultra-light ALPs explored in Ref. [32]. In fact these two mixing scenarios, in-source and 
in the IGMF, are complementary to each other and cover a huge range of ALP mass. De- 
tection of source signatures can be used to constrain the ALP mass as well as to put limit 
on the IGMF. Indeed, there are hints from recent studies of ultra high-energy cosmic ray 
data and TeV blazars that the IGMF can be much smaller than a nG [66], in which case 
the polarization and spectral signatures of in-source photon-ALP mixing that we explored 



will not be destroyed. Photon-ALP mixing in the ~ /uG galactic magnetic field over kpc 
coherence length scale is also negligible. 

5. Conclusions 

Axions and axion-like particles appear in many extentions of the standard model of particle 
physics. Photon-axion/ALP mixing in the presence of an external electromagnetic field 
constitutes one of the most exploited signals for astrophysical and laboratory axion and 
ALP searches. Gamma-ray bursts are the most powerful source of keV-MeV photons 
in nature, which are most probably synchrotron radiation. These photons originate and 
propagate inside the GRB jet with high magnetic field. We have shown that strong photon- 
ALP conversion takes place in GRB jet in the ~ 100 keV observed energy range, distorting 
the standard synchrotron polarization pattern. We have also shown that when the magnetic 
field direction in the photon propagation coherence length is perpendicular to the magnetic 
field direction in the synchrotron radiating region, the photon energy spectrum will be 
steeper than the expected spectrum from synchrotron theory, thus providing an explanation 
for the anomalous spectra of ~ 20% of the observed GRBs. We found that the photon- ALP 
conversion occurs within a large range of possible GRB and ALP parameters, being almost 
independent of the ALP mass for sufficiently small ALP masses (m a < 10 -6 eV). Further 
modification due to mixing in the intergalactic magnetic field is not expected in case the 
IGMF is < 1 nG and/or the ALP mass is > 10~ 14 eV. 

Large statistics expected to be collected by a number of future missions that are 
devoted to measure GRB polarization in the keV-MeV range will be crucial to search for 
ALP signals due to their mixing with photons inside GRBs. 
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